Digital twinning of the human ventricular activation sequence to Clinical 12-lead ECGs and magnetic resonance imaging using realistic Purkinje networks for in silico clinical trials.

Cardiac in silico clinical trials can virtually assess the safety and efficacy of therapies using human-based modelling and simulation. These technologies can provide mechanistic explanations for clinically observed pathological behaviour. Designing virtual cohorts for in silico trials requires exploiting clinical data to capture the physiological variability in the human population. The clinical characterisation of ventricular activation and the Purkinje network is challenging, especially non-invasively. Our study aims to present a novel digital twinning pipeline that can efficiently generate and integrate Purkinje networks into human multiscale biventricular models based on subject-specific clinical 12-lead electrocardiogram and magnetic resonance recordings. Essential novel features of the pipeline are the human-based Purkinje network generation method, personalisation considering ECG R wave progression as well as QRS morphology, and translation from reduced-order Eikonal models to equivalent biophysically-detailed monodomain ones. We demonstrate ECG simulations in line with clinical data with clinical image-based multiscale models with Purkinje in four control subjects and two hypertrophic cardiomyopathy patients (simulated and clinical QRS complexes with Pearson's correlation coefficients > 0.7). Our methods also considered possible differences in the density of Purkinje myocardial junctions in the Eikonal-based inference as regional conduction velocities. These differences translated into regional coupling effects between Purkinje and myocardial models in the monodomain formulation. In summary, we demonstrate a digital twin pipeline enabling simulations yielding clinically consistent ECGs with clinical CMR image-based biventricular multiscale models, including personalised Purkinje in healthy and cardiac disease conditions.

Cardiac Corpuscles: A "New" Morphofunctional Entity? / Corpúsculos Cardíacos: ¿Una "Nueva" Entidad Morfofuncional?

SUMMARY: The existence of "transitional muscular structures" between subendocardial branches (Purkinje fibers) and ventricular working muscle fibers (WF) was first described by the German anatomist, Kurt Goerttler, in 1964. He designated them as "subendocardial nucleus organs." He supposed such fibers functioned as mechanoreceptors, controlling of the intensity of contraction of the ventricular musculature. Brazilian anatomist Ferraz de Carvalho described similar structures in 1993. A thorough literature search failed to identify any other research articles confirming or denying their existence. The objective of this work was to find such structures in subendocardial ventricular walls in human hearts. We collected fifteen formalin-preserved hearts from the Anatomy Department of São Paulo University and sectioned the apical portions on the right and left ventricles according to method used by Goerttler. We utilized conventional histology (light microscopy- LM), scanning electron microscopy (SEM), and a new preservation method called micro- plastination (MP). At the anterior wall of the right ventricle in the subendocardial region between the interventricular septum and moderator band, we found several bundles of fusiform and helicoidal fibers of similar histology to the WF. The bundles measured between 400 and 1150 µm in length and were separated from adjacent muscular fibers by thin collagen fiber, thus acting as a "pseudo capsule." Some structures seemed to be linked to PF and were appeared to be lymphatic and blood vessels and nerves. We called those structures "cardiac corpuscles" (CC). The observation of the previously "unknown" CC in this initial study confirmed the previous descriptions and its discovery may contribute to new perspectives in the study of cardiac muscle structure and function.

La existencia de "estructuras musculares de transición" entre los ramos subendocárdicos (fibras de Purkinje) y las fibras musculares ventriculares activas(FMV) fue descrita por primera vez por el anatomista alemán Kurt Goerttler en 1964, quien las denominó "órganos del núcleo subendocárdico". Supuso que tales fibras funcionaban como mecanoreceptores, controlando la intensidad de la contracción de la musculatura ventricular. El anatomista brasileño Ferraz de Carvalho describió estructuras similares en 1993. Una búsqueda bibliográfica exhaustiva no logró identificar ningún otro artículo de investigación que confirmara o negara su existencia. El objetivo de este trabajo fue encontrar dichas estructuras en las paredes ventriculares subendocárdicas de corazones humanos. Recolectamos 15 corazones conservados en formalina del Departamento de Anatomía de la Universidad de São Paulo y seccionamos las porciones apicales de los ventrículos derecho e izquierdo según el método utilizado por Goerttler. Utilizamos histología convencional (microscopía de luz-LM), microscopía electrónica de barrido (SEM) y un nuevo método de conservación llamado microplastinación (MP). En la pared anterior del ventrículo derecho en la región subendocárdica entre el tabique interventricular y la banda moderadora, encontramos varios haces de fibras fusiformes y helicoidales de histología similar a la FMV. Los haces medían entre 400 y 1150 µm de longitud y estaban separados de las fibras musculares adyacentes por una fina fibra de colágeno, actuando así como una "pseudocápsula". Algunas estructuras parecían estar vinculadas a la fibras de purkinje y parecían ser vasos linfáticos, sanguíneos y nerviosos. Llamamos a esas estructuras "corpúsculos cardíacos" (CC). La observación del CC previamente "desconocido" en este estudio inicial confirmó las descripciones anteriores y su descubrimiento puede contribuir a nuevas perspectivas en el estudio de la estructura y función del músculo cardíaco.

Histological study on the heart ventricle of Egyptian bovines (Bos aegyptiacus).

Background: The heart ventricles have thicker walls than atrium as they pump blood through blood vessels into all body organs. Aim: This study aimed to describe the histological changes of the heart ventricles in Egyptian bovine (Bos aegyptiacus) with special reference to Purkinje fibers. Methods: A total of 10 male Egyptian bovines of 1-10 years old were divided into three groups according to age; immature, mature, and adult animals. Results: The histological sections from all examined animals' groups revealed three different layers of the wall of both right and left ventricles; endocardium, myocardium, and epicardium. The endocardium was lined with endothelium and filled with fibrous connective tissue. The endocardium of adult bovine was the thickest. Purkinje fibers appeared of pale cytoplasm with few myofibrils. They were present in the deep layer of the endocardium and in the myocardium. The size of Purkinje fibers and the amount of their myofibrils appeared to be increased with advanced age. Bundles of cardiac muscles were the main constituent of the myocardium. The myocardial bundles were separated by fine connective tissue in immature animals that showed an increased amount in the adult animals. The hypereosinophilic cardiac muscle cells were observed in the ventricles of both mature and adult animals suggesting hypercontraction during rigor mortis. An external layer of the ventricles was the epicardium which consisted of connective tissue and covered with mesothelium. Conclusion: Overall, this study revealed histological changes in the wall of the ventricle and Purkinje fibers of Egyptian bovines (B. aegyptiacus) in relation to age. Additionally, the hypereosinophilia of the cardiac muscle cells was recorded in the ventricles of mature and adult bovines.

Purkinje physiology and pathophysiology.

There has always been an appreciation of the role of Purkinje fibers in the fast conduction of the normal cardiac impulse. Here, we briefly update our knowledge of this important set of cardiac cells. We discuss the anatomy of a Purkinje fiber strand, the importance of longitudinal conduction within a strand, circus movement within a strand, conduction, and excitability properties of Purkinjes. At the cell level, we discuss the important components of the ion channel makeup in the nonremodeled Purkinjes of healthy hearts. Finally, we discuss the role of the Purkinjes in forming the heritable arrhythmogenic substrates such as long QT, heritable conduction slowing, CPVT, sQT, and Brugada syndromes.

Crista Supraventricularis Purkinje Network and Its Relation to Intraseptal Purkinje Network.

Using transparent specimens with a dual color injection, microscopy, and computer tomography, this report shows that the right and left ventricular subendocardial Purkinje networks are connected by an extensive septal network in the bovine heart. The septal network is present along the entire septum except at a free zone below ventricular valves. Being the only communication of the basal right septum with the right free wall, the supraventricular crest is an enigmatic but not, by any means, hidden muscular structure. It is one of the last structures to be activated in human heart. It is shown here that the supraventricular crest Purkinje network connects the anterosuperior right ventricular basal free wall Purkinje network to anterior right ventricular basal septal Purkinje network. It is suggested that the stimulus initiated at middle left ventricular endocardium will activate the supraventricular crest. The intraseptal connection found between the basal left ventricular subendocardial septal Purkinje network and the right ventricular basal septal Purkinje network is, probably, the pathway for the stimulus. An anatomic basis is provided to explain why the inflow tract contracts earlier than the outflow tract in the right ventricle systole. Anat Rec, 2017. © 2017 Wiley Periodicals, Inc. Anat Rec, 300:1793-1801, 2017. © 2017 Wiley Periodicals, Inc.

Morphological study of the atrioventricular conduction system and Purkinje fibers in yak.

We studied the morphology of the atrioventricular conduction system (AVCS) and Purkinje fibers of the yak. Light and transmission electron microscopy were used to study the histological features of AVCS. The distributional characteristics of the His-bundle, the left bundle branch (LBB), right bundle branch (RBB), and Purkinje fiber network of yak hearts were examined using gross dissection, ink injection, and ABS casting. The results showed that the atrioventricular node (AVN) of yak located in the right side of interatrial septum and had a flattened ovoid shape. The AVN of yak is composed of the slender, interweaving cells formed almost entirely of the transitional cells (T-cells). The His-bundle extended from the AVN, and split into left LBB and RBB at the crest of the interventricular septum. The LBB descended along the left side of interventricular septum. At approximately the upper 1/3 of the interventricular septum, the LBB typically divided into three branches. The RBB ran under the endocardium of the right side of interventricular septum, and extended to the base of septal papillary muscle, passed into the moderator band, crossed the right ventricular cavity to reach the base of anterior papillary muscle, and divided into four fascicles under the subendocardial layer. The Purkinje fibers in the ventricle formed a complex spatial network. The distributional and cellular component characteristics of the AVCS and Purkinje fibers ensured normal cardiac function.

Cellular Physiology and Clinical Manifestations of Fascicular Arrhythmias in Normal Hearts.

Fascicular ventricular arrhythmias represent a spectrum of ventricular tachycardias dependent on the specialized conduction system. Although they are more common in structurally abnormal hearts, there is an increasing body of literature describing their role in normal hearts. In this review, the authors present data from both basic and clinical research that explore the current understanding of idiopathic fascicular ventricular arrhythmias. Evaluation of the cellular electrophysiology of the Purkinje cells shows clear evidence of enhanced automaticity and triggered activity as potential mechanisms of arrhythmias. Perhaps more importantly, heterogeneity in conduction system velocity and refractoriness of the left ventricular conduction system in animal models are in line with clinical descriptions of re-entrant fascicular arrhythmias in humans. Further advances in our understanding of the conduction system will help bridge the current gap between basic science and clinical fascicular arrhythmias.

Protective Role of False Tendon in Subjects with Left Bundle Branch Block: A Virtual Population Study.

False tendons (FTs) are fibrous or fibromuscular bands that can be found in both the normal and abnormal human heart in various anatomical forms depending on their attachment points, tissue types, and geometrical properties. While FTs are widely considered to affect the function of the heart, their specific roles remain largely unclear and unexplored. In this paper, we present an in silico study of the ventricular activation time of the human heart in the presence of FTs. This study presents the first computational model of the human heart that includes a FT, Purkinje network, and papillary muscles. Based on this model, we perform simulations to investigate the effect of different types of FTs on hearts with the electrical conduction abnormality of a left bundle branch block (LBBB). We employ a virtual population of 70 human hearts derived from a statistical atlas, and run a total of 560 simulations to assess ventricular activation time with different FT configurations. The obtained results indicate that, in the presence of a LBBB, the FT reduces the total activation time that is abnormally augmented due to a branch block, to such an extent that surgical implant of cardiac resynchronisation devices might not be recommended by international guidelines. Specifically, the simulation results show that FTs reduce the QRS duration at least 10 ms in 80% of hearts, and up to 45 ms for FTs connecting to the ventricular free wall, suggesting a significant reduction of cardiovascular mortality risk. In further simulation studies we show the reduction in the QRS duration is more sensitive to the shape of the heart then the size of the heart or the exact location of the FT. Finally, the model suggests that FTs may contribute to reducing the activation time difference between the left and right ventricles from 12 ms to 4 ms. We conclude that FTs may provide an alternative conduction pathway that compensates for the propagation delay caused by the LBBB. Further investigation is needed to quantify the clinical impact of FTs on cardiovascular mortality risk.

Incidence, Distribution and Morphology of Left Ventricular False Tendons in Cat Hearts.

The incidence, distribution, and macro- and microscopic structures of left ventricular false tendons (LVFTs) in the cat heart were studied using 25 normal and 57 diseased hearts. The fibrous bands were observed in the left ventricle of all 82 cat hearts examined and most commonly extended between the papillary muscles and the ventricular septum. Histologically, the LVFTs were composed of central Purkinje fibres and surrounding dense collagenous fibres covered by endothelium. There was no appreciable difference in the incidence, distribution or morphology of LVFTs between the normal and the diseased hearts, indicating that LVFTs are a common anatomic variant in the cat heart.

Estimation of Purkinje trees from electro-anatomical mapping of the left ventricle using minimal cost geodesics.

The electrical activation of the heart is a complex physiological process that is essential for the understanding of several cardiac dysfunctions, such as ventricular tachycardia (VT). Nowadays, patient-specific activation times on ventricular chambers can be estimated from electro-anatomical maps, providing crucial information to clinicians for guiding cardiac radio-frequency ablation treatment. However, some relevant electrical pathways such as those of the Purkinje system are very difficult to interpret from these maps due to sparsity of data and the limited spatial resolution of the system. We present here a novel method to estimate these fast electrical pathways from the local activations maps (LATs) obtained from electro-anatomical maps. The location of Purkinje-myocardial junctions (PMJs) is estimated considering them as critical points of a distance map defined by the activation maps, and then minimal cost geodesic paths are computed on the ventricular surface between the detected junctions. Experiments to validate the proposed method have been carried out in simplified and realistic simulated data, showing good performance on recovering the main characteristics of simulated Purkinje networks (e.g. PMJs). A feasibility study with real cases of fascicular VT was also performed, showing promising results.

Ungulates heart model: a study of the Purkinje network using India ink injection, transparent specimens and computer tomography.

The Purkinje network is not macroscopically visible in human hearts. Sunao Tawara found himself in trouble in the early 1900s, when studying the human heart network. He gained a much better understanding of the net after starting to work with ungulates' hearts. The ungulate heart is proposed as an auxiliary didactic model for the study of the human conduction system. This work provides a detailed description of the India ink injection technique to allow a naked eye visualization of the Purkinje network. The heart muscle was made diaphanous for direct visualization of the ungulate heart intramyocardial network, and computer tomography was employed for visualization of the three dimensional structure of the whole network. The intramyocardial network in the interventricular septum was identified. The pattern of the Purkinje network is described as a connected noneulerian graph, and its possible implications on the mechanism of arrhythmias is discussed. The main differences between the ungulate and human heart conduction systems are stressed.

Computational generation of the Purkinje network driven by clinical measurements: the case of pathological propagations.

To properly describe the electrical activity of the left ventricle, it is necessary to model the Purkinje fibers, responsible for the fast and coordinate ventricular activation, and their interaction with the muscular propagation. The aim of this work is to propose a methodology for the generation of a patient-specific Purkinje network driven by clinical measurements of the activation times related to pathological propagations. In this case, one needs to consider a strongly coupled problem between the network and the muscle, where the feedback from the latter to the former cannot be neglected as in a normal propagation. We apply the proposed strategy to data acquired on three subjects, one of them suffering from muscular conduction problems owing to a scar and the other two with a muscular pre-excitation syndrome (Wolff-Parkinson-White). To assess the accuracy of the proposed method, we compare the results obtained by using the patient-specific Purkinje network generated by our strategy with the ones obtained by using a non-patient-specific network. The results show that the mean absolute errors in the activation time is reduced for all the cases, highlighting the importance of including a patient-specific Purkinje network in computational models.

Functional, anatomical, and molecular investigation of the cardiac conduction system and arrhythmogenic atrioventricular ring tissue in the rat heart.

BACKGROUND: The cardiac conduction system consists of the sinus node, nodal extensions, atrioventricular (AV) node, penetrating bundle, bundle branches, and Purkinje fibers. Node-like AV ring tissue also exists at the AV junctions, and the right and left rings unite at the retroaortic node. The study aims were to (1) construct a 3-dimensional anatomical model of the AV rings and retroaortic node, (2) map electrical activation in the right ring and study its action potential characteristics, and (3) examine gene expression in the right ring and retroaortic node. METHODS AND RESULTS: Three-dimensional reconstruction (based on magnetic resonance imaging, histology, and immunohistochemistry) showed the extent and organization of the specialized tissues (eg, how the AV rings form the right and left nodal extensions into the AV node). Multiextracellular electrode array and microelectrode mapping of isolated right ring preparations revealed robust spontaneous activity with characteristic diastolic depolarization. Using laser microdissection gene expression measured at the mRNA level (using quantitative PCR) and protein level (using immunohistochemistry and Western blotting) showed that the right ring and retroaortic node, like the sinus node and AV node but, unlike ventricular muscle, had statistically significant higher expression of key transcription factors (including Tbx3, Msx2, and Id2) and ion channels (including HCN4, Cav3.1, Cav3.2, Kv1.5, SK1, Kir3.1, and Kir3.4) and lower expression of other key ion channels (Nav1.5 and Kir2.1). CONCLUSIONS: The AV rings and retroaortic node possess gene expression profiles similar to that of the AV node. Ion channel expression and electrophysiological recordings show the AV rings could act as ectopic pacemakers and a source of atrial tachycardia.

Characterization and modeling of the peripheral cardiac conduction system.

The development of biophysical models of the heart has the potential to get insights in the patho-physiology of the heart, which requires to accurately modeling anatomy and function. The electrical activation sequence of the ventricles depends strongly on the cardiac conduction system (CCS). Its morphology and function cannot be observed in vivo, and therefore data available come from histological studies. We present a review on data available of the peripheral CCS including new experiments. In order to build a realistic model of the CCS we designed a procedure to extract morphological characteristics of the CCS from stained calf tissue samples. A CCS model personalized with our measurements has been built using L-systems. The effect of key unknown parameters of the model in the electrical activation of the left ventricle has been analyzed. The CCS models generated share the main characteristics of observed stained Purkinje networks. The timing of the simulated electrical activation sequences were in the physiological range for CCS models that included enough density of PMJs. These results show that this approach is a potential methodology for collecting knowledge-domain data and build improved CCS models of the heart automatically.

Anatomical and molecular mapping of the left and right ventricular His-Purkinje conduction networks.

Functioning of the cardiac conduction system depends critically on its structure and its complement of ion channels. Therefore, the aim of this study was to document both the structure and ion channel expression of the left and right ventricular His-Purkinje networks, as we have previously done for the sinoatrial and atrioventricular nodes. A three-dimensional (3D) anatomical computer model of the His-Purkinje network of the rabbit heart was constructed after staining the network by immunoenzyme labelling of a marker protein, middle neurofilament. The bundle of His is a ribbon-like structure and the architecture of the His-Purkinje network differs between the left and right ventricles. The 3D model is able to explain the breakthrough points of the action potential on the ventricular epicardium during sinus rhythm. Using quantitative PCR, the expression levels of the major ion channels were measured in the free running left and right Purkinje fibres of the rabbit heart. Expression of ion channels differs from that of the working myocardium and can explain the specialised electrical activity of the Purkinje fibres as suggested by computer simulations; the expression profile of the left Purkinje fibres is more specialised than that of the right Purkinje fibres. The structure and ion channel expression of the Purkinje fibres are highly specialised and tailored to the functioning of the system. The His-Purkinje network in the left ventricle is more developed than that in the right ventricle and this may explain its greater clinical importance.

Inter-model consistency and complementarity: learning from ex-vivo imaging and electrophysiological data towards an integrated understanding of cardiac physiology.

Computational models of the heart at various scales and levels of complexity have been independently developed, parameterised and validated using a wide range of experimental data for over four decades. However, despite remarkable progress, the lack of coordinated efforts to compare and combine these computational models has limited their impact on the numerous open questions in cardiac physiology. To address this issue, a comprehensive dataset has previously been made available to the community that contains the cardiac anatomy and fibre orientations from magnetic resonance imaging as well as epicardial transmembrane potentials from optical mapping measured on a perfused ex-vivo porcine heart. This data was used to develop and customize four models of cardiac electrophysiology with different level of details, including a personalized fast conduction Purkinje system, a maximum a posteriori estimation of the 3D distribution of transmembrane potential, the personalization of a simplified reaction-diffusion model, and a detailed biophysical model with generic conduction parameters. This study proposes the integration of these four models into a single modelling and simulation pipeline, after analyzing their common features and discrepancies. The proposed integrated pipeline demonstrates an increase prediction power of depolarization isochrones in different pacing conditions.

Fast multiscale modeling of cardiac electrophysiology including Purkinje system.

In this paper, we present a modeling methodology to couple the cardiac conduction system to cardiac myocytes through a model of Purkinje-ventricular junctions to yield fast and realistic electrical activation of the ventricles. A patient-specific biventricular geometry is obtained from processing computed tomography scan data. A one-manifold implementation of the fast marching method based on Eikonal-type equations is used for modeling heart electrophysiology, which facilitates the multiscale 1-D-3-D coupling at very low computational costs. The method is illustrated in in-silico experiments where we analyze and compare alternative pacing strategies on the same patient-specific anatomy. We also show very good agreement between the results from the proposed approach and more detailed and comprehensive biophysical models for modeling cardiac electrophysiology. The effect of atrioventricular delay on the distribution of activation time in myocardium is studied with two experiments. Given the reasonable computational times and realistic activation sequences provided by our method, it can have an important clinical impact on the selection of optimal implantation sites of pacing leads or placement of ablation catheter's tip in the context of cardiac rhythm management therapies.

Estimation of local orientations in fibrous structures with applications to the Purkinje system.

The extraction of the cardiac Purkinje system (PS) from intensity images is a critical step toward the development of realistic structural models of the heart. Such models are important for uncovering the mechanisms of cardiac disease and improving its treatment and prevention. Unfortunately, the manual extraction of the PS is a challenging and error-prone task due to the presence of image noise and numerous fiber junctions. To deal with these challenges, we propose a framework that estimates local fiber orientations with high accuracy and reconstructs the fibers via tracking. Our key contribution is the development of a descriptor for estimating the orientation distribution function (ODF), a spherical function encoding the local geometry of the fibers at a point of interest. The fiber/branch orientations are identified as the modes of the ODFs via spherical clustering and guide the extraction of the fiber centerlines. Experiments on synthetic data evaluate the sensitivity of our approach to image noise, width of the fiber, and choice of the mode detection strategy, and show its superior performance compared to those of the existing descriptors. Experiments on the free-running PS in an MR image also demonstrate the accuracy of our method in reconstructing such sparse fibrous structures.

Flexible modeling for anatomically-based cardiac conduction system construction.

We present a method to automatically deploy the peripheral section of the cardiac conduction system in ventricles. The method encodes anatomical information thorough rules that ensure that Purkinje network structures generated are realistic and comparable to those observed in ex-vivo studies. The core methodology is based in non-deterministic production rules that are parameterized by means of statistical functions. Input parameters allow the construction of a great diversity of Purkinje structures that could be incorporated in fine element ventricular models to perform electrophysiology simulations. Resulting Purkinje trees show good geometrical approximations of Purkinje core network and bundles when compared to histological diagrams and do not require user interaction. Simulations carried out with these models result in activation sequences remarkably similar to micro-electrode electrical mapping studies.

The anatomy of the cardiac conduction system.

All the myocytes within the heart have the capacity to conduct the cardiac impulse. A population of myocytes is specialized so as to generate the cardiac impulse and then to conduct it from the atrial to the ventricular chambers. This population has become known as the conduction system. Anatomists who seek to demonstrate the location of the components of this system must contend with the fact that the components of the system cannot be distinguished from the working myocardial elements by gross dissection. In important presentations to the German Pathological Society in 1910, rules were suggested for the histological distinction of these conducting cells. These rules proposed that the myocytes, to be considered as part of the conduction system, should be histologically discrete, traceable from section to section in serially prepared material, and if to be considered as tracts, should be insulated by fibrous tissue from the adjacent myocytes. Immunohistochemical techniques have now been developed that better demonstrate the distinction between the cells specialized to conduct from working myocytes. These new techniques, for the most part, confirm the accuracy of the initial descriptions. They also reveal additional areas with the characteristics of conduction tissues. These additional areas are located in a paranodal area adjacent to the sinus node, in the vestibules of both atrioventricular valvar orifices, and in a partial ring around the aortic root. In this review, we describe all these features, emphasizing the relationship of the newly recognized components to the established parts of the cardiac conduction system, and how the new findings need to be assessed in the light of the old criteria.